Resonantly enhanced pair production in a simple diatomic model 
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CsJ '■ Abstract 



A new mechanism for the production of electron-positron pairs from the interaction of a laser 
field and a fully stripped diatomic molecule in the tunneling regime is presented. When the laser 



^ ■ field is turned off, the Dirac operator has resonances in both the positive and the negative energy 



continua while bound states are in the mass gap. When this system is immersed in a strong laser 



^ ■ field, the resonances move in the complex energy plane: the negative energy resonances are pushed 

o: 

■ to higher energies while the bound states are Stark shifted. It is argued here that there is a pair 
, production enhancement at the crossing of resonances by looking at a simple 1-D model: the nuclei 

od 

■ are modeled simply by Dirac delta potential wells while the laser field is assumed to be static and 
(N" 

of finite spatial extent. The average rate for the number of electron-positron pairs produced is 
evaluated and the results are compared to the single nucleus and to the free cases. It is shown that 
5-H I positrons are produced by the Resonantly Enhanced Pair Production (REPP) mechanism, which 

is analogous to the resonantly enhanced ionization of molecular physics. This phenomenon could 
be used to increase the number of pairs produced at low field strength, allowing the study of the 
Dirac vacuum. 

PACS numbers: 31.50.Gh,12.20.Ds,03.75.Dg,02.30.Gp 
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There has been a tremendous amount of efforts in the last few decades to increase laser 
intensities and it is now conceivable to reach laser intensities above 10^^ W/cm^ such 
that the corresponding field strength is comparable to Coulomb fields in matter. For slightly 
lower intensities, this has led to the development of nonperturbative models of ionization such 
as tunnelling ionization in atomic physics and Charge Resonance Enhanced Ionization 
(CREI) in molecular physics j^. The new regime of laser intensities now available allows 
the study of new physical effects where relativistic and Quantum Electrodynamics (QED) 
corrections start to be important Q]. Among the QED effects, one of the most important 



phenomenon is the long-sought Schwinger's mechanism js], which consists of the decay of 
the vacuum of a static electric field into electron-positron pairs. In the Dirac interpretation, 
this can be seen as a tunneling of electrons from the negative energy sea to the positive 
continuum. This has never been observed experimentally because it requires field strength 
on the order of Es ~ 10^® V/m, which are not available experimentally (it corresponds to 
a laser intensity of 10^^ W/cm^). However, given the new experimental advances and the 
novel laser technologies, there has been a renewed interest in this process and new ideas 
have emerged which could allow to probe the QED vacuum. Therefore, many variations 



Jl5|. A 



16|. 



of Schwinger's original idea using different field configurations were proposed 
semi-classical theory of relativistic tunneling ionization has also been considered 

In this letter, a new mechanism is proposed to enhance pair production from a laser 
field interacting with a fully stripped diatomic molecule, which we denote the Resonantly 
Enhanced Pair Production (REPP) process. This mechanism is analogous to the CREI 
process, which is well-known in molecular physics [3|. It proceeds in the following way and 
is depicted in Fig. [H First, let us consider the case when the electric field F is zero and look 
at the spectrum of the Dirac operator for a system having two nuclei. Also, for simplicity, we 
consider only two bound states: the ground state and the first excited state. If the nuclear 
charge is not too large {Z < 137), the bound states are in the mass gap, in the energy range 
[—1710^,1710^]. On the other hand, the negative and positive energy continua incorporate 
the so-called Ramsauer-Townsend resonances (RTR) (these are shown in dark on Fig. [T]), 
related to the backscattering of waves on the potential wells. This results in a peak-valley 
structure in the spectral density function in both continua When the electric field is 
turned on and the interatomic distance R is varied, the position of resonances in the complex 
energy plane changes: the RTR of the negative energy sea are pushed to higher energies 



while the RTR of the positive energy continuum are pulled towards lower energies [17|] . The 
bound states become resonances (they gain an imaginary part and thus, becomes unstable 
states) and are Stark-shifted by A ~ -^FR/2. When a RTR of the negative energy sea 
crosses a resonance from the bound states or a RTR from the positive energy continuum, 
the transition from the negative to the positive energy continuum is enhanced, yielding 
a higher pair production rate (channel 1). The same mechanism occurs for the excited 
state: when it crosses a RTR from the positive energy sea, the ionization of the molecule 
is enhanced, facilitating also the transition from the negative energy sea to the resonance 
because it reduces the Pauli blocking (if the electron is ionized, the excited state resonance 
is "free" and can receive a new electron that tunneled from the negative energy state). This 
is the channel 2. Therefore, if the field or the interatomic distance is not too large, the pair 
production occurs by a two-steps process: electrons from the negative energy sea tunnel to 
one of the bound state resonances, which is then followed by an ionization process where the 
same electron tunnels to the positive energy continuum. This produces a flux of electron 
and positrons, propagating in opposite direction. There is another possibility which occurs 
when the RTR from the positive and negative continua cross each other (channel 3). Then, 
it is possible to have a direct transition from the negative to the positive energy continua, 
which also enhances pair production but which effect is usually more modest than other 
channels. These effects are important at large atomic distance. For small atomic distances, 
such as in those achieved in Heavy ion collisions, the wave functions of each nucleus overlap 
and another mechanism is responsible for pair production: the effective electric charge in 
this case approaches 2Z in which case the ground state has a lower energy level. In the 
electric field, it is then easier for an electron to tunnel from the negative energy sea and 
again, this can lead to enhanced pair production. 

To confirm these ideas, a calculation of the pair production rate is performed in a very 
simple 1-D model: the nuclei are modeled by Dirac potential wells while the laser field 
is considered in the adiabatic limit and is taken as a static electric field. In a previous 
3-D model of the one electron molecule in superintense laser fields, it was shown how 



to adapt Volkov solutions of the Dirac equation to the two centre problem [18|. In the 
present work, we use a simpler model which allows us to understand the mechanisms and 
the physics behind the pair production from the interaction of lasers with fully stripped 
diatomic molecules. We show that REPP can be utilized to enhance pair production at low 
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FIG. 1: Description of the REPP mechanism. The dark grey regions are the position of 
resonances (high density of states) while the hght grey regions correspond to accessible 

energies which have a lower density of state. 



field strength and large internuclear distance. This conclusion is reached by conducting a 
comparative study with the cases where there is no nuclei (Schwinger's process) and where 
there is only one nucleus. 

b calculate the average rate of pair produced, we follow the discussion presented in 



19 



- |21| and assume that the electric field vanishes at x = ±00 and thus, that it has a finite 
extent in space. In this case, it is possible to define the "asymptotic states" at x = =Foo: in 
these regions, the particles are free and there is a natural separation between the negative 
and positive energy solutions. This allows us to evaluate the number of pairs produced from 
a solution of the Dirac equation. It should be noted here that boundary conditions on the 
wave function at x = =Foo are obtained from the time- dependent case where the Lehmann- 
Symanzik-Zimmermann (LSZ) asymptotic conditions, which implies the vanishing of the 



field at t = =Foo, can be used 22|]. In the time-independent case however, the latter cannot 
be fulfilled directly: we have to consider localized wave packets which are effectively in the 
free region when t = =Foo. From these considerations, it is possible to evaluate the average 




FIG. 2: Simple model to study REPP. The electric field has a finite extent in space. In 
blue is the Klein region where it is possible to have a transition from a negative to positive 

energy states. 



number of pairs (n) and it has been argued in 19N21| that this observable is given by 



where A is the coefficient of the positive energy solution propagating towards x = +00, at 
the right of the potential (see Fig. [2]). This formula is valid for a time- independent external 
field where solutions are labeled by energy. Similar formula have been derived in {23, 24 1. 

Thus, the calculation of pair production reduces to a transmission-refiection problem 
where the incident, refiected and transmitted waves are given respectively by: 

^inc.(a;)=^;(p)e^^(^)^ (2) 
^ref.(a;)=5^^(-p)e-^^(^)^ (3) 
Ar.nsXx) = Au{k)e'''^'^^'. (4) 

Here, u, v are the positive and negative energy free spinors: their explicit expression will be 
given below. 

We are considering the 1-D Dirac equation, which is given by 

E'iIj{x) = [-icaA + o^mc^ + Aq{x) + V{x)\ ^(a;), (5) 
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where E is the eigenenergy, c is the hght velocity, m is the electron mass and iplx) is the two- 
components spinor wave function. The electric potential is divided in three spatial regions 
as (see Fig. [2]) 



Ao(x) 



2FL for X e (-00, -L] 

-F{x-L) for a; G(-L,L) 
for X G [L, 00) 



(6) 



where F is the field strength (we are working in a gauge where the vector potential is = 
and the electric field is related to the potential as E^ = —dxAo{x)) and 2L is the length over 
which the electric field is constant. Outside the interval [—L,L], the electric field vanishes. 
For the other scalar potential, three cases are considered: 

1. No nucleus: V{x) =0 



2. Single nucleus: V{x) 



-g6{x) 



3. Two nuclei: V{x) = —gS{x — i?) — g6{x + R) 



where g is the potential strength (physically, it is related to the charge of the nucleus). It 



was shown in 17| that the Dirac delta potential wells can be characterized by the following 
boundary conditions: \im^^oip{( + e) = ^i^e^oG'ipiC " ^) (C is the potential well position), 
which relates the wave function on the right and the left of the potential well. The matrix 
components are given by G12 = G21 = and 



G 



1 



11 



1 + ^ 



9 , .9 
4c2 c 



G. 



22- 



(7) 



Thus, the wave function between the potential wells is a solution of Eq. ([5]) without V{x). 
Then, Eqs. ([7]) is used to match the wave function at x = 0, ±R, for case 2 and 3 respectively. 
The Dirac equation can be solved analytically by decoupling the two spinor components 



and by letting y{x) = e y ^ y ^ j ■ Then, the Dirac equation becomes a system of 
equations with well-known solutions in terms of parabolic cylinder functions f/(7, z) 25| : 



'ljj{x) = CiUa{x) + C2Ub{x), 
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(8) 



where ci,2 are integration constants and where we defined 

Ua,lix) = f/(7,y(x)), 

Un2ix) = mc 



^ e''-^U{^ + l,y{x)), 



Ub,2{x) 



mc 



2F 
2F 



-e '4 



U{-j-l,y{x)). 



(9) 
(10) 

(11) 
(12) 



Here, we have 7 



i. The last missing ingredient to compute pair production are 



2F 2- 

the negative and positive energy free spinors. There is a well-known procedure to compute 
these quantities where one seek plane wave solutions. This yields 



u{k) 



v{p) 



'2E 



+ ck{E) 



1 



^E - ck{E) 

^{E - 2FL) + cp{E) 



^2{E - 2FL) 



■y/{E-2FL)-cp{E) 



(13) 
(14) 



where k{E) = \^/E'^ - m?c^ and p{E) = \^ {E ~ 2FLY - m?c^. 

Now, imposing the continuity of the wave function at the region boundaries and using 
Eq. ([7]) to include the potential wells, we obtain the following conditions for the three cases: 

1. No nucleus: 

v{p)e''P'^ + Bv{-p)e-'^'' = a^Ua{~L) + a^U^i-L), 
aiUa{L) + a2Ub{L) = Au{k)e'^'' . 



2. Single nucleus: 



;(p)e'^'^ + Bv{-p)e-'^'' = a^Ua{-L) + a^U^i-L), 



aiUaiO) + a2f/fe(0) = [hUaiO) + hUi 



ipL 



3. Two nuclei: 



hUaiL) + b2UhiL) = Auik)e 



v{p)e'P'' + Bv{-p)e-'P'' = a,Ua{~L) + a^U^i-L), 
aiUa{-R) + a^Uki-R) = G'' [hUa{-R) + h^-R)] , 

hiUa{R) + hUb{R) = [c^Ua{R) + C2Ub{R)] , 

Cit/,(L) + C2Ub{L) = Au{k)e''P'^. 



These are systems of equations which can be used to solve for the integration constants 
A, B, ai 2, &i,2, Ci,2- The numerical results are presented in the following. 

The particle spectrum d{n) / dEdt as a function of interatomic distance is plotted in Fig. 
[3] for the two nuclei case and for g = 0.8 (the ground state energy of this potential well 
corresponds to the energy of the Is orbital of the U^^"*" ion). The other parameters are 
chosen as L = 100., F = 0.09, in units where h = c = m = 1 and e = y/a. In this figure, 
there is a clear enhancement of pair production at the position of resonances. Moreover, as 
R is varied, it is possible to see how the resonances are moving in the energy plane: some 
of the RTR in the negative energy sea are moving to higher energies while the ground state 
and the first excited state are Stark shifted. More interesting is the fact that when the RER 
crosses with the bound states, there is an enhancement of pair production (channel 1 or 2). 

In Fig. m the total rate d{n)/dt as a function of the interatomic distance is presented and 
compared to the other cases. This figures show clearly that the number of pairs is enhanced 
by the REPP at larger R: the position of peak in the rate corresponds to the interatomic 
distance where the ground state resonance crosses with RTR (channel 1). Thus, the number 
of pairs is approximately an order of magnitude higher than other cases and therefore, 
REPP is an important process for pair production in laser-matter interaction. The largest 
enhancement however occurs at small R where it reaches two order of magnitude above the 
single nucleus case. This suggests that an experiment using HIC in a laser field could also 
be used to probe the Dirac vacuum. 
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